Global phase synchronization in an array of time-delay systems 
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We report the identification of global phase synchronization (GPS) in a linear array of unidi- 
rectionally coupled Mackey-Glass time-delay systems exhibiting highly non-phase-coherent chaotic 
attractors with complex topological structure. In particular, we show that the dynamical organi- 
zation of all the coupled time-delay systems in the array to form GPS is achieved by sequential 
synchronization as a function of the coupling strength. Further, the asynchronous ones in the ar- 
ray with respect to the main sequentially synchronized cluster organize themselves to form clusters 
before they achieve synchronization with the main cluster. We have confirmed these results by 
estimating instantaneous phases including phase difference, average phase, average frequency, fre- 
quency ratio and their differences from suitably transformed phase coherent attractors after using a 
nonlinear transformation of the original non-phase-coherent attractors. The results are further cor- 
roborated using two other independent approaches based on recurrence analysis and the concept of 
localized sets from the original non-phase-coherent attractors directly without explicitly introducing 
the measure of phase. 

PACS numbers: 05.45.Xt,05.45.Pq 



I. INTRODUCTION 



Chaotic phase synchronization (CPS), referred to as 
the locking of the phases of the coupled chaotically evolv- 
ing dynamical systems, has been investigated in ensem- 
bles of globally coupled arrays [l|-0, [IlTtl3| . networks 
of oscillators [IJ-QjI wit h ap plications to electrochem- 
istry |7|, [8| , laser systems fig ) , 1 3 1 , cardiorespirator y sys- 



iMg], ecology HJig, cli- 



tems [18H20II. neuroscience 

matology [27H29J ] . etc. Even though the notion of CPS 
is well studied in low dimensional systems, there ex- 
ist very little indepth studies in higher dimensional sys- 
tems such as time-delay systems which are essentially 
infinite-dimensional in nature and often exhibit high- 
dimensional, highly non-phase-coherent hyperchaotic at- 
tractors with complex topological structure. Conse- 
quently, estimating phase explicitly to identify phase syn- 
chronization in such systems is quite difficult. Neverthe- 
less, CPS has been recently demonstrated in two cou- 
pled piecewise linear and Mackey-Glass time-delay sys- 
tems [3fj, |3l| by introducing a nonlinear transformation 
of the original dynamical variable to recast the original 
non-phase-coherent hyperchaotic attractors into smeared 
limit cycle-like attractors in order to facilitate the esti- 
mation of the phase variable using the available meth- 
ods. However, these investigations are carried out so 
far only in two coupled time-delay systems. In view 
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of the widespread applications of CPS in ensembles of 
coupled oscillators, here we investigate the existence of 
global phase synchronization (GPS) in an array of unidi- 
rectionally coupled Mackey-Glass time-delay systems and 
analyse the mechanism behind the dynamical organiza- 
tion of the coupled oscillators to form GPS. At first, we 
use the nonlinear transformation introduced in [30L |3l| 
to estimate explicitly the phases of all the oscillators in 
the array and identify the existence of GPS in the ar- 
ray. Further, the existence of GPS is also confirmed 
from the original non-phase-coherent chaotic attractors 
themselves using two independent approaches, namely 
recurrence analysis [32|, [33| and the concept of localized 
sets [11]. 

In addition, we will show that the onset of GPS in 
such arrays does not happen instantaneously, but instead 
takes place as a form of sequential synchronization. For 
lower values of coupling strengths the phases of nearby 
systems get already synchronized with the drive system 
in contrast to the far away systems. This sequential syn- 
chronization of chaotic systems can have applications in 
communication systems |35| . Furthermore, other non- 
synchronized time-delay systems with respect to the se- 
quentially synchronized cluster display clusters of phase 
synchronized states among themselves before they be- 
come synchronized with the large cluster in the sequence 
to form global phase synchronization. This clustering is 
observed when the group of oscillators splits into sub- 
groups such that all the oscillators within one cluster 
move in perfect phase synchrony. This clustering is con- 
sidered to be particularly significant in biological sys- 
tems [36l438| . Recently cluster synchronization in an ar- 
ray of three chaotic lasers without delay was reported [39| 
as well. Also global synchronization via cluster formation 



1.5 
1.2 
0.9 
0.6 
0.3 




2 
1.5 

1 
0.5 






0.3 0.6 0.9 1.2 1.5 0.3 0.6 0.9 1.2 




0.3 0.6 0.9 1.2 1.5 0.3 0.6 0.9 1.2 

x 15 (t+%) x lg (t+x) 



FIG. 1: (a) Non phase-coherent chaotic attractor of the drive 
system given by Eq. (QJ,). (b) Transformed attractor (using 
Eq. ((2|) of the drive system, (c-f) Transformed attractors of 
some randomly selected response systems (i — 8, 12, 15, 19) in 
the projected phase space (xi(t + r),Zi(t + r)), where they 
look like smeared limit cycle attractors in the absence of the 
coupling along with the Poincare points represented by filled 
circles. Here r has been chosen as r = 20.0. 



has been observed in coupled phase oscillators without 
time-delay [10f with simultaneous phase slips of all oscil- 
lators, where quantized phase shifts in these phase slips 
have been observed. By increasing the coupling, a bi- 
furcation tree from high-dimensional quasiperiodicity to 
chaos to quasiperiodicity and periodicity has also been 
found. 

The remaining paper is organized as follows: In Sec. ITII 
we will describe briefly the coupling configuration and 
the nature of chaotic attractors exhibited by the Mackey- 
Glass time-delay system. The existence of global phase 
synchronization from the transformed attractors is dis- 
cussed in Sec. IIII[ which is also confirmed from the orig- 
inal non-phase-coherent chaotic attractors using recur- 
rence analysis and the concept of localized sets in Sec. lIVI 
Finally, we summarize our results in Sec. [Vj 



II. LINEAR ARRAY OF MACKEY-GLASS 
TIME-DELAY SYSTEMS 

The Mackey- Glass time-delay system was originally de- 
duced as a model for blood production in patients with 
leukemia |4fJ, and it has been well studied in the lit- 




FIG. 2: The first four maximal Lyapunov exponents X ma x of 
the Mackey-Glass time-delay system fl^) for the parameter 
values q = 0.1,/? = 0.2, r £ (14,37) in the absence of the 
coupling C. 



erature for its hyperchaotic behavior [4lH43| and has 
also been exp erimentally realized using analog electronic 
circuits [4J|. In this paper, we consider a linear array 
of unidirectionally coupled Mackey-Glass systems with 
free-end boundary conditions represented by the follow- 
ing system of coupled nonlinear first order ordinary dif- 
ferential equations, 



ii{t) 
±i(t) 
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a>iXi(t 



{l + x^t-rYY 
a % Xi(t - t) 



(la) 



{1 + Xi(t 
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t = 2,3,.-.,JV,(lb) 



where, a, j3, c are the system parameters, r is the time- 
delay and C is the coupling strength. We have fixed 
the parameter values at a.\ = 0.2, f3 = 0.1, c = 10.0, 
t = 20.0 and the values of the nonlinear parameter a.% of 
the response systems in the array are chosen randomly 
in the range a, € (0.17,0.20), so that all the systems 
are effectively nonidentical. For our simulations, we have 
fixed the total number of oscillators in the array as N — 
20, though we confirmed our results for N — 50 also (see 
Appendix A). 

The uncoupled drive (pj,) exhibits a highly non-phase- 
coherent chaotic attractor for the chosen parameter val- 
ues, which is depicted in Fig. [1^. The first four largest 
Lyapunov exponents of the uncoupled drive system are 
shown in Fig. H as a function of the delay time r. 
Note that the phase calculated directly from the origi- 
nal non-phase-coherent chaotic attractor (Fig. QJ,) can- 
not yield monotonically increasing behaviour as it has 
several closed loops, which also contribute to the phase 
information, other than the main center of rotation of the 
major part of the trajectories [30, |31(. To overcome this 
problem a nonlinear transformation is introduced so as 
to rescale the original non-phase-coherent chaotic attrac- 
tor into smeared limit cycle-like attractor with a single 
center of rotation. The transformation is performed by 
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FIG. 3: (a-d) Phase differences (A(/>i,i — 4>i — 4>i) of the randomly selected systems (i = 8,12,15,19) from the array of coupled 
Mackey-Glass time-delay systems for different values of the coupling strength C. A more detailed description can be found in 
the text. 



introducing a new state variable [30t |31|. 

z(t + r) = x(t)x(t + t)/x{t + t), 



(2) 



where f is the optimal value of time-delay to be cho- 
sen in order to avoid any additional center of rotation. 
This functional form of transformation (along with a de- 
lay time f) has been identified by generalizing the trans- 
formation used in the case of chaotic attractors in the 
Lorenz systems [l|. Now, the projected trajectory in the 
new state space (x(t + r),z(t + r)) (Fig. [1}d) resembles 
that of a smeared limit cycle-like attractor with a sin- 
gle fixed center of rotation. It is also to be noted that, 
even though the transformed attractor has sharp turns 
in the vicinity of the common center, it does not have 
any closed loops as in the original non-phase-coherent 
attractor. Otherwise, the transformed attractor would 
not give rise to monotonically increasing phase resulting 
in exact matching of the phases of the coupled time-delay 
systems |30l[3l1|. 



III. GPS FROM THE TRANSFORMED 
ATTRACTOR 

In this section, we will show that the global phase syn- 
chronization in the array of Mackey-Glass time-delay sys- 
tems ([T]) is attained by a sequential phase synchronization 
of the oscillators in the array as the coupling strength is 



increased. Further we will also demonstrate that the re- 
maining non-synchronized oscillators in the array form 
synchronized clusters among themselves before attaining 
GPS. 

We use the same nonlinear transformation ([2]) to recast 
the original non-phase-coherent chaotic attractors of all 
the N = 20 oscillators into smeared limit cycle-like at- 
tractors. We have fixed the value of the optimal value f 
in Eq. ^ as f = 8.0 for all the N oscillators. Instead one 
can also choose different values for f for different oscilla- 
tors, as they are nonidentical systems with a parameter 
mismatch, to obtain more exact unfolding for different 
attractors. However, we find f — 8.0 for all the oscilla- 
tors is adequate for our purpose in the following study. 
We have calculated the instantaneous phases of all the os- 
cillators using the Poincare section technique [J, [2| from 
their corresponding transformed attractors. Projected 
trajectories of randomly selected response systems (i — 
8,12,15,19) in the array (JTJd) into the new state space 
(xi(t +r),Zi(t + t)), where they look like smeared limit 
cycle-like attractors with a fixed center of rotation, are 
shown in Figs. [TJc-f). Filled circles in these plots corre- 
spond to the Poincare points. 

Phase differences, A(f>i = (f>i — (f>i, between the drive 
and some randomly selected response systems (i — 
8, 12, 15, 19 ) in the array (JTJd) are shown in Figs. [3] for 
different values of the coupling strength. They increase 
monotonically in the absence of coupling (C = 0.0) in- 
dicating that all the oscillators are in an asynchronous 
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FIG. 4: The average phase difference (r](t)) for different values 
of the coupling strength C. For C = 0.0 the phases of all 
the systems are unbounded so the phase difference increases 
linearly with time but for C — 1.2 the phases of all the systems 
are bounded, showing a high quality phase synchronization. 



state. Phase slips in the phase differences for small val- 
ues of the coupling strength indicates that the oscillators 
are in the transition state to GPS. Further increase in the 
value of the coupling strength results in a strong bound- 
edness of the phases of the oscillators. For sufficiently 
large C, the phase differences become zero (Figs. [3]) in- 
dicating the existence of phase synchronization between 
the drive and the response systems. It is evident from 
the Figs. [3] that the 8th oscillator is synchronized with 
the drive at C = 0.7, while the other systems are in the 
transition state, whereas 12th oscillator is synchronized 
with drive only at C = 0.9. The other two oscillators 
with the index i — 15 and i = 19 reach synchronization 
with the drive for further larger values of the coupling 
strength, C — 1.0 and 1.2, respectively. Therefore, it is 
clear that the nearby oscillators to the drive system in 
the array are synchronized first as the coupling strength 
is increased implying that the global phase synchroniza- 
tion is reached by sequential phase synchronization of the 
coupled oscillators in the array. To confirm the existence 
of GPS, we have calculated the average phase difference, 
r](t), defined as 



v(t) = 



N 



1 N 



(3) 



j=2 



The average phase difference (r](t)) for different values 
C is shown in Fig. Has a function of time t. In the 
absence of the coupling (C = 0.0), the phases of all the 
oscillators evolve independently and hence their average 
phase difference increases linearly with time. Further 
increase in C induces the entrainment of phases of the 
oscillators and at the value of coupling strength C = 
1.2 the average phase difference of all the N oscillators 
becomes exactly zero, showing a high quality GPS in the 
array. 

The emergence of GPS through a sequential phase syn- 
chronization is also characterized by calculating the time 
averaged phase ({4>i)) and the time averaged frequency 
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FIG. 5: (a) Time averaged phase ({4>i}) and (b) time averaged 
frequency ((fii)) of all the systems plotted as a function of the 
system index i for various values of the coupling strength C. 



((Oj)) of each of the oscillators in the array which are 
defined as 



(Mt)) = (2wk + 2n f \ ) , (4<i*<4 + 0a) 



(Qi(t)) = lim - / <t>i(t)dt 



T^ooT 



(4b) 



where t\ is the time of the k crossing of the flow with 
the Poincare section of the i th attractor and (...)t denotes 
a time average. The average phase and the average fre- 
quency are shown in Figs. [5ji and Figs. \Sjp, respectively, 
for different values of the coupling strength as a func- 
tion of the oscillator index. A random distribution of 
the average phase (Fig. EK(i)) and the average frequency 
(Fig.[S}3(i)) for the value of the coupling strength C = 0.1 
indicate that the coupled oscillators in the array evolve 
almost independently. A slight increase in the coupling 
strength (to C — 0.3) results in synchronous evolution of 
the first 5 oscillators in the array as seen in Fig. [5ji(ii) 
and Fig.[5p(ii). For C = 0.5, Fig. rji(ih) and Fig.[5p(iii) 
indicate that the first 6 oscillators are synchronized. It is 
also to be noted from these figures that the other desyn- 
chronized oscillators form synchronized clusters among 
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FIG. 6: (Colour online) Snap shots of the node vs node diagrams (that is oscillator index vs oscillator index plots) indicating 
the sequential phase synchronization and the organization of cluster states for different values of coupling strength. The 
different symbols indicate that the corresponding nodes are phase synchronized, (a) non phase synchronized case for C = 0.1 
(b) First four oscillators in Eq. ([lb) are phase synchronized with the drive system for C = 0.27. (c), (d) and (e) Sequential 
phase synchronization and the formation of small cluster states for C = 0.5, 0.82 and 1.08, respectively, and (f) Global phase 
synchronization for C — 1.2. 



themselves. In particular, the oscillators with the indices 
8—11 synchronize among themselves to form a cluster, 
while the oscillators with the indices 13 — 14, and 19 — 20 
form separate small clusters. These clusters can also be 
clearly visualized by plotting the oscillator index plots 
as we will illustrate below. It is also to be noted that 
even when the total number of oscillators in the array 
is increased, the phenomenon remains qualitatively the 
same though the sizes of the clusters will increase appro- 
priately (see Appendix A below for some details for N = 
50). We also note that the results remain qualitatively 
unaltered even for different sets of random values for the 
nonlinear parameters, (%, confirming the robustness of 
our results. A similar transition to PS through clustering, 
termed as hard transition for large coupling strength, in a 
chain of diffusively coupled Rossler oscillators with large 
frequency mismatch have been observed [9j but in the 
periodic state due to the suppression of chaotic attrac- 
tors. For further larger values of the coupling strength, 
the desynchronized oscillators form similar small clusters 
among themselves before attaining GPS. The average 
phase and the average frequency illustrated in Fig.[5^,(iv) 
and Fig. [5]3(iv) for C — IX) indicates that most of the 
nearest oscillators are synchronized with the drive, while 
the oscillators with the indices 18 — 19 form a small sep- 
arate cluster. All the oscillators in the array become 
phase/frequency locked and evolve in synchrony (GPS) 
with each other for the coupling strength C = 1.12 as de- 
picted in Fig.[S^,(v) and Fig. [S}d(v) and they continue to 
be in a stable phase/frequency synchronized state which 



is shown in Fig.[5^(vi) and Fig. [Sb(vi) for C = 2.0. 

The mechanism for the formation of clusters and GPS 
may be explained as follows. Due to the mismatches in 
the nonlinear parameters, o^, all the individual systems 
in the array evolve independently with different phases 
(phase mismatches), and correspondingly with frequency 
mismatches, for small values of the coupling strength C. 
As C is increased further, the oscillators with nearest 
frequencies in the array synchronize first to form clus- 
ters among themselves leaving the clusters with relatively 
large frequency mismatch in isolation (see Figs.ltjlandrT3"|). 
Further increase in C results in the formation of a single 
large cluster whose constituents exhibit a coherent phase 
oscillation with the drive due to the decomposition of the 
clusters away from the drive in the array. Consequently 
GPS results in the system. Similar mechanism has been 
identified in ensemble of coupled Rossler oscillators with 
frequency mismatches [45j (without time-delay). 

The above dynamical organization of GPS via sequen- 
tial phase synchronization and the clustering can be also 
be visualized clearly by using snap shots of the oscillators 
in the index vs index plot, as node vs node diagrams, as 
shown in Fig. [6] The oscillators that evolve with identical 
values of the average phase/frequency are assigned with 
identical shapes. The diagonal line in Fig. [5^ for C = 0.1 
corresponds to the oscillator index i = j and they evolve 
independently. Figure |6Jd indicates that the first four os- 
cillators in the array are synchronized with the drive for 
C = 0.27. As discussed above, three small clusters are 
seen in Fig. [BJ: for C = 0.5 while the first 6 oscillators 
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FIG. 7: (Colour online) (a) Time averaged phase ({&}) and 
(b) Time averaged frequency ({f2j)), i = 1,2, ••• ,20 plotted 
as a function of the coupling strength C £ (0, 1.5). Here for 
each value of C we have plotted the average phase/frequency 
of all the N = 20 oscillators which is shown by the filled circles. 
Insets show some of the systems get synchronized themselves 
to form small clusters (each subgroup of oscillator is differ- 
entiated by different types of circles) before they synchronize 
with the drive system to form GPS. 
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FIG. 8: (Colour online) (a) The frequency difference 
(Afiij, j = 2, 3, • • • ,N) and (b) the frequency ratio 
(Qj/ili, j — 2, 3, • ■ • , N) are plotted as a function of the 
coupling strength C £ (0, 1.5). Each line corresponds to the 
difference/ratio between a response system and the drive sys- 
tem. The black filled circles indicate the average frequency 
difference/ratio of all the (N — 1) response systems from the 
drive system. 



form a large synchronized cluster. Similar small clusters 
are shown in Figs. [Hi and [His for C = 0.82 and 1.08, re- 
spectively, in addition to the single large cluster formed 
by sequential phase synchronization of the oscillators in 
the array. Finally, GPS of all the oscillators in the array 
is illustrated in Fig.[5J for C = 1.2. 

For a global picture of the emergence of GPS, we have 
plotted the average phase ((0i)) and the average fre- 
quency ((fij)) of all the TV oscillators as a function of 
the coupling strength C in Figs. [7J There is an absence 
of any correlation among the average phases (Fig.[?k.) and 
the average frequencies (Fig. [7b) of different oscillators 
for low values of the coupling strength as revealed by the 
random distributions of their values. The random distri- 
bution of (<pi) and (Oj) are organized into few clusters 
for C > 0.5 as may be evident from the Figs. [7] Global 
phase synchronization emerges for C > 1.12 as may be 
seen from the insets. Small synchronized clusters formed 
by the remaining asynchronous oscillators for larger val- 
ues of C can also be appreciated from the insets. 

We have also plotted the frequency difference 
(Afiij, j = 2, 3, •••,JV) and the frequency ratio 
(ftj/Cli, j = 2, 3, • • • , N) of all the oscillators with that 
of the drive as a function of the coupling strength by dif- 
ferent types of lines in Fig. [Sb- and Fig. [Bb, respectively. 



The black filled circles in Fig. [8^ and Fig. [8b correspond 
to the average frequency difference and the average fre- 
quency ratio of all the oscillators with that of the drive. 
The substantial saturation in their values for C > 1.12 
indicates the emergence of GPS. 

The emergence of global phase synchronization in the 
array can also be quantified using the well-known order 
parameter 47] . 



Re 



itfj 



1 N 



(5) 



where 4>j(t) denotes the instantaneous phase of the j th 
system, ip(t) is the average phase and (...)* denotes a time 
average. If all the systems are in a phase synchronized 
state then R sw 1. The order parameter (R) is plotted 
in Fig. [HI for the number of oscillators N — 20 and N 
— 50 as a function of the coupling strength C. As C is 
increased, R also increases and for N = 20 the critical 
value of the coupling is C > 1.12 and for TV — 50, C will 
be > 2.4, for the value of R w 1 confirms the existence 
of GPS in the array of coupled time-delay systems. 



IV. GPS FROM THE ORIGINAL 
NON-PHASE-COHERENT ATTRACTOR 

In this section, we demonstrate the existence of GPS 
from the original non-phase-coherent chaotic attractors 
using two different app roaches, namely recurrence quan- 
tification analysis [32, [33| and the concept of localized 
sets [34| without estimating explicitly the measure of 
phase. 



A. GPS using recurrence analysis 

Several measures of complexity which quantify small 
scale structures in the recurrence plots have been pro- 
posed and are known as recurrence quantification anal- 
ysis (RQA) [32|, |33|. Certain measures have also been 
introduced to characterize and identify different kinds of 
synchronization transitions in coupled chaotic systems. 
These measures have the advantage of applicability in 
the analysis of experimental systems and, in particular, 
in the case of very small available data sets. Further, 
these measures can also be used in the case of non-phase- 
coherent chaotic/hyperchaotic attractors of time-delay 
systems 30, 31, 46], where it is difficult and often even 
impossible to calculate the phase explicitly. Among the 
available recurrence quantification measures, we use the 
Correlation of Probability of Recurrence (CPR) and the 
generalized autocorrelation function P{t) to confirm the 
existence of GPS in the array of coupled time-delay sys- 
tems ([J), both qualitatively and quantitatively. 

A criterion to quantify phase synchronization between 
two systems is the Correlation of Probability of Recur- 
rence (CPR) defined as 



CPR=(Pi(t)P 2 (t))/ai*2, 



(6) 



where P(t) is the generalized autocorrelation function 
represented as 



P(t) 



1 



N-t 



N -t 



^eie-wXi-Xi+tW), 



(7) 



where O is the Heaviside function, Xi is the i data 
point of the system X, e is a predefined threshold, ||.|| is 
the Euclidean norm, and N is the number of data points, 
Pi, 2 means that the mean value has been subtracted and 
ai t 2 are the standard deviations of Pi{t) and -P2W, re- 
spectively Looking at the coincidence of the positions of 
the maxima of P (t) of the systems, one can qualitatively 
identify PS [32|, [33|. If both systems are in CPS, the 
probability of recurrence is maximal at the same time 
t and CPR re 1. If they are not in CPS, the maxima 
do not occur simultaneously and hence one can expect a 
drift in both the probability of recurrences resulting in 
low values of CPR. 

The generalized autocorrelation function of the drive 
P\{t) and that of some response systems (i — 8, 12, and 
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FIG. 9: The phase order parameter (R) for the number of 
oscillators N — 20 and N = 50 as a function of the coupling 
strength indicating global phase synchronization in the array 
of coupled Mackey-Glass time-delay systems. 



19) Pg(t), Pi 2 (t), and P 19 (t) are depicted in Fig. [TO] for 
different values of the coupling strength. In the absence 
of coupling (C = 0.0), all systems evolve independently 
and hence the maxima of their respective generalized au- 
tocorrelation functions do not occur simultaneously as 
shown in Fig. 110a . On increasing the coupling strength, 
oscillators with a lower value of index in the array become 
synchronized first resulting in sequential phase synchro- 
nization and this can also be identified from the gener- 
alized autocorrelation functions of the response systems 
in the array. For instance, Ps{t), Pii{t), and Pig(i) are 
shown along with P\{t) in Fig. [TOb for C = 0.4. It is 
clear from this figure that the maxima of the drive P\ (t) 
and those of the response Ps(t) are in complete agree- 
ment with each other (Fig. [TObi) indicating the existence 
of PS between them. On the other hand, only some of 
the maxima of the response system P\i (t) are in coinci- 
dence with those of the drive fFig.lTObii) illustrating that 
the response system i = 12 is in transition to PS, whereas 
the maxima of the response system P19 (t) do not coincide 
with those of the drive (Fig. [TObiii') indicating that the 
response system i = 19 is in an asynchronous state for the 
same value of C . For C — 1.2, almost all of the positions 
of the peaks of the generalized auto correlation functions 
Pi(t), Ps(t), Pi2(t), and Pig(t) are in agreement with 
each other as illustrated in Fig. [TOT c) confirming the ex- 
istence of GPS via sequential phase synchronization. It 
is also to be noted that the magnitudes of the peaks of 
all the oscillators have generally of different values and 
the differences in the heights of the peaks indicate that 
there is no correlation in the amplitudes of the coupled 
systems. Furthermore, the formation of clusters by the 
other asynchronous oscillators in the array can also be 
realized by plotting their respective generalized autocor- 
relation functions, which will show that all their maxima 
are in good agreement with each other, whereas there 
exists a drift between them and the maxima of the se- 
quentially synchronized cluster. 

The existence of GPS via sequential phase synchroniza- 
tion is also quantified using value of the index CPR of the 
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FIG. 10: (Colour online) Generalized autocorrelation functions of the drive Pi(t) and randomly selected response systems 
(i = 8, 12, and, 19) Ps{t), Pi2{t), and Pig(i) indicating (a) Non-phase-synchronization for C — 0.0, (b) generalized autocorrelation 
functions for C = 0.4 (bi) PS between the systems 1 and 8, (bii) approximate PS between the systems 1 and 12 and (biii) non 
PS between the systems 1 and 19, and (c) PS between all the systems (i = 1,8, 12, and 19) for C — 1.2. 



response systems with the drive as shown in Fig. 1111 The 
different lines correspond to the index of the oscillators 
(i = 2,8,12,19) in the array. It is evident from the fig- 
ure that the oscillators with increasing index attain the 
value of unity in a sequence as a function of the coupling 
strength and finally for C > 1.12 the CPR of all the re- 
sponse systems with the drive reaches unity confirming 
that all the coupled oscillators are in GPS. The mean 
value of CPR of all the response systems in the array is 
shown as filled circles, which also confirms the existence 
of GPS for C> 1.12. 
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B. GPS using the concept of localized sets 

Recently, an interesting framework to identify CPS, 
namely, the concept of localized sets 34| has been in- 
troduced. This approach provides an easy and efficient 
way to detect CPS especially in complicated non-phase- 
coherent attractors. The basic idea of this concept is to 
define a typical event in one of the systems and then ob- 
serve the other system whenever this event occurs. These 
observations give raise to a set D. Depending upon the 
property of this set D, one can state whether PS exists 
or not. The coupled systems evolve independently if the 
sets obtained by observing the corresponding events in 
the systems spread over the attractor of the systems. On 



FIG. 11: (Colour online) The index CPR as a function of the 
coupling strength C. Different lines correspond to the CPR 
of different (i — 2, 8, 12, and 19) response systems with the 
drive system. The filled circles correspond to the mean value 
of the CPR of all the (N — 1) systems in the array. 



the other hand, if the sets are localized on the attractors 
then CPS exists between them. 

We have confirmed the existence of the GPS in the 
linear array of Mackey- Glass time-delay systems (dJ by 
using this concept of localized sets. Now, we will demon- 
strate the existence of GPS via sequential phase synchro- 
nization in the randomly selected response systems (i = 
1,8,12,19). We have defined the event as Poincafe sec- 



C = 0.0 



C = 0.4 



C = 0.6 



C = 0.76 



C = 1.2 



i = l 



i= 12 



19 



1.4 

1.2 

-P 1 

+ 

i 0.8 

'"' 0.6 

0.4 



1.4 
1.2 
1 
0.8 
0.6 
0.4 



1.4 
1.2 

? > 
^ 0.8 

* 0.6 

0.4 



1.4 
1.2 
1 
0.8 
0.6 
0.4 



•".- 4 'V, 4 ., 



(a) 



1.4 | 

1.2 
1 
0.8 
0.6 
0.4 (e) 



0.4 0.6 0.8 1 1.2 1.4 
x,(t) 






(b) -'^ 

0.4 0.6 0.8 1 1.2 1.4 
xjt) 




0.4 0.6 0.8 1 1.2 1.4 
x 12 (t) 



ft - ?. 






. ; •'.' 


't, 


'- !«.. 


** ,•% 


•it 




(d) 







0.4 0.6 0.8 1 1.2 1.4 




1.4 
1.2 

P 1 

+ 

i 0.8 

"' 0.6 

0.4 




1.4 

1.2 

P 1 

X 0.8 

"' 0.6 

0.4 



0.4 0.6 0.8 1 1.2 1.4 
x,(t) 



0.4 0.6 0.8 1 1.2 1.4 

Xj(t) 



0.4 0.6 0.8 1 1.2 1.4 
x,d<) 



0.4 0.6 0.8 1 1.2 1.4 
x lg (t) 




On) 




-^ 



1.4 
1.2 

P 1 

+ 

i 0.8 

H 0.6 

0.4 



0.4 0.6 0.8 1 1.2 1.4 
x,(t) 



1.4 

1.2 

P 1 

4 0.8 

* 0.6 

0.4 



0.4 0.6 0.8 1 1.2 1.4 




0.4 0.6 0.8 1 1.2 1.4 
x 12 (t) 



1.4 

1.2 

P 1 

4 0.8 

*" 0.6 

0.4 



0.4 0.6 0.8 1 1.2 1.4 
x is (t) 




(q) 



0.4 0.6 0.8 1 1.2 1.4 
x,(t) 




« 

0.4 0.6 0.8 1 1.2 1.4 





1.4 




1.2 


+ 

K 1 


1 
0.8 
0.6 




0.4 




0.4 0.6 0.8 1 1.2 1.4 
X 12 (t) 




0.4 0.6 0.8 1 1.2 1.4 

x Ig (t) 



FIG. 12: (Colour online) First row (a-q) corresponds to the attractors of the drive system (i = 1) and the rows (b-r, c-s, and 
d-t) correspond to the attractors of some randomly selected response systems (i = 8, 12, 19). The ' + ' marks represent the 
events (Poincare sections) in the corresponding attractors. In (a-d) the sets (represented by the filled circles) are spread over 
the attractors and hence there is no CPS for the value of coupling strength C — 0.0. In (e-h) for C — 0.4 and in (i-1, m-p, q-t) 
the sets are localized confirming the existence of GPS in the array for C = 0.6, 0.76, and 1.2, respectively. 



tions in the attractors indicated as ' + ' marks in Figs. 1121 
The set, indicated as filled circles, obtained by observing 
the drive system (i — 1) whenever the defined event oc- 
curs in the response system (i ~ 8) is shown in Fig. [12] 
(a) and that obtained by observing the response systems 
i = 8,12,19 whenever the defined event occurs in the 
drive system are shown in Figs. Q2] (b-d) for the value 
of coupling strength C — 0.0. As the obtained sets are 
spread over the attractors, all the systems evolve inde- 
pendently and there is no CPS in the absence of coupling 
between them. Further when we increase the coupling 
strength to C = 0.4, the oscillator (i = 8) is partially 
synchronized with the drive as the sets are almost local- 
ized but the sets in the oscillators i = 12,19 are spread 
over the attractors which means that they are not yet 
phase synchronized with the drive system. This is shown 
in Figs[T2](e-h). Again increasing the coupling strength 
to C = 0.6, the sets are further bounded to a small region 
over the attractors which shows that the oscillators i — 
8,12 are synchronized with the drive, but the oscillator i 
— 19 is less phase synchronized with the drive which is 
represented by the spread of the events over the attractor 



as shown in Figs[T2](i-l)- Further, the FigsfT2l (m-p) and 
FigsfT2](q-t) indicate the situation for C = 0.76 and C — 
1.2, respectively, where all the oscillators are now phase 
synchronized with the drive as the sets are localized over 
the attractor confirming the existence of GPS in an ar- 
ray via sequential phase synchronization as the coupling 
strength is increased. 

Further, the formation of clusters can also be realized 
using the concept of localized sets by defining the event 
among the response systems that form the clusters and 
observing the other response systems that are in the same 
cluster. In this case the obtained sets by observing the 
event in the drive will spread over the attractor of the 
response systems, while the sets obtained by observing 
the event among the response systems that form a cluster 
will be localized on their respective attractors. 



V. SUMMARY AND CONCLUSION 

We have demonstrated the existence of global phase 
synchronization via sequential phase synchronization 
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FIG. 13: (Colour online) Snap shots of the node vs node diagram indicating the sequential phase synchronization and the 
organization of cluster states of TV = 50 oscillators for different values of coupling strength, C. The different symbols indicate 
that the corresponding nodes are phase synchronized, (a) Non-phase synchronized case for C = 0.1, (b) First seven oscillators 
in Eq. (JTb) are phase synchronized with the drive system for C = 0.4, (c), (d) and (e) Sequential phase synchronization and 
the formation of small cluster states for C = 0.53, 1.9 and 2.3, respectively, and (f)Global phase synchronization of N = 50 
oscillators for C = 2.5. 



in an array of coupled Mackey- Glass time-delay sys- 
tems with parameter mismatches, which exhibit highly 
non-phase-coherent attractors with complex topological 
structure. Further, we have also shown that the re- 
maining asynchronous systems will organize themselves 
to form different clusters before they get phase synchro- 
nized with the main cluster to form global phase synchro- 
nization. We have confirmed the existence of GPS via se- 
quential phase synchronization by estimating the phase 
difference as a function of the coupling strength, the aver- 
age frequency and the average phase as a function of the 
oscillator index and the coupling strength after calculat- 
ing the phase variables from the transformed attractors. 
We have also demonstrated the existence of sequential 
phase synchronization and the formation of clusters by 
the remaining oscillators using the average frequency, av- 
erage phase and specifically using the index vs index plot 
of the oscillators. Furthermore, we have also confirmed 
the existence of GPS via sequential phase synchroniza- 
tion using the recurrence quantification measures and 
the concept of localized sets which are calculated from 
the original non-phase-coherent attractors of the coupled 
Mackey-Glass time-delay systems. It is also to noted that 
we have obtained similar transitions to GPS via cluster- 
ing even in the hyperchaotic regimes of Fig. 2 for the 
Mackey-Glass systems and also in the coupled piecewise 
linear time-delay systems (with five positive Lyapunov 



exponents) . 



Acknowledgments 

The work of R.S. and M.L. has been supported by the 
Department of Science and Technology (DST), Govern- 
ment of India sponsored IRHPA research project, and 
DST Ramanna program of M.L. D.V.S. has been sup- 
ported by the Alexander von Humboldt Foundation. J.K. 
acknowledges the support from EU under project No. 
240763 PHOCUS(FP7-ICT-2009-C). 



Appendix A: 

The dynamical organization of GPS via sequential 
phase synchronization and the cluster formation can be 
visualized by using the snap shots of N = 50 oscillators 
in the index vs index plot as shown in Fig|13l The diag- 
onal line in FigJT3] (a) for the value of coupling C = 0.1 
corresponds to the oscillator index i = j and the oscilla- 
tors evolve independently. Further in Fig[l3] (b) , the first 
seven oscillators in the array are synchronized with the 
drive and the oscillators 8 — 11 form a separate cluster for 
the coupling strength C = 0.4. Further, the first eleven 



11 



oscillators in the array are synchronized and four small 
separate clusters are seen in FigJT3](c) for C = 0.53. Sim- 
ilar small clusters are formed in Fig[l3] (d) for C = 1.9 
while the first thirty one oscillators form a large synchro- 
nized cluster. Further, in Fig[13] (e) the first forty four 



oscillators are synchronized with the drive and the oscil- 
lators 45 — 50 form a separate cluster for C — 2.3. Finally, 
the occurrence of GPS of all the oscillators in the array is 
illustrated in Figfnj] (f ) for the value of coupling strength 
C = 2.5. 
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